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Abstract. It is controversial what is the true role of entanglement in two-photon 
virtual-state spectroscopy [Saleh et al, Phys. Rev. Lett. 80, 3483, 1998], a two-photon 
absorption spectroscopic technique that can retrieve information about the energy level 
structure of an atom or a molecule. The consideration of closely related techniques, 
such as multidimensional pump-probe spectroscopy [Roslyak et al, Phys. Rev. A 
79, 063409, 2009], might suggest that spectroscopic information retrieved in the two- 
photon absorption process is the same regardless of the classical or quantum nature of 
the light source. Here, we solve this debate by making use of a full quantum formalism 
to show that the ability to obtain information about the energy level structure of 
a medium requires the existence of temporal (frequency) correlations between the 
absorbed photons. Moreover, we show that these correlations are not the only requisite 
for retrieving such information. In fact, it is a combination of both, the presence of 
frequency correlations and its specific spectral shape, which makes the realization 
of two-photon virtual-state spectroscopy possible. This result helps clarifying the 
discussion whether entanglement is needed or not, and also, to specify the type of 
two-photon source that needs to be used in order to experimentally perform the two- 
photon virtual-state spectroscopy technique. 
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1. Introduction 

The process of two-photon absorption (TPA), the light-induced transition between two 
energy levels of a medium mediated by the absorption of two photons, is a building 
block of some technologies aimed at probing the structure of atoms and molecules, such 
as two-photon microscopy [1] and two-photon spectroscopy [2] . In particular, nonlinear 
two-photon spectroscopy has become an invaluable tool [3] , where the capability of TPA 
is exploited to obtain information about a sample that would not be accessible otherwise. 

With the advent of light sources that generate entangled photon pairs [3], new 
phenomena in TPA processes have been unveiled. The linear dependence of the 
two-photon absorption rate on the photon flux [5], two-photon induced transparency 
[6], virtual-state spectroscopy [7J [8] and the selectivity of double-exciton states of 
chromophore aggregates [9] are effects that have been attributed to the presence of 
entanglement. However, the link between entanglement and the new effect observed 
is sometimes blur. So, might not be the ultimate cause of some of these effects an 
accompanying characteristic unrelated to its entangled nature? This is the case, for 
instance, of the cancelation of dispersion effects which can take place in optical coherence 
tomography [TOlfTT]. the observation of remote dispersion cancelation [12] and temporal 
modulation [13] of entangled photons, effects that could also be produced by frequency- 
correlated photons, which nonetheless might be non-entangled [T3] , 

Here, we consider one important spectroscopic application whose capabilities have 
been associated to the use of entangled photons pairs, namely two-photon virtual-state 
spectroscopy. The absorption of two photons by an atom or a molecule induces a 
transition between two of its energy levels that match the overall energy of the incident 
photons. The quantum mechanical calculation of the TPA transition probability shows 
that its value can be understood as a weighted sum of many energy non- conserving 
atomic transitions (virtual-state transitions) [T5l 116] between energy levels. 

The virtual-state transitions, a signature of the medium, can be revealed 
experimentally by introducing a delay between the two absorbed photons, and averaging 
over different experimental realizations with different temporal correlations between the 
photons [7]. Can we retrieve the sought-after information (energy level structure) with 
any type of frequency correlations between the photons? It has been suggested that 
spectroscopic information resident in the TPA signal in multidimensional pump-probe 
spectroscopy [XT] is essentially the same, regardless of the classical or quantum nature 
of the light used. As stated recently in [9j, it remains however an open question, to 
what extent these effects constitute genuine entanglement effects and whether they can 
be reproduced, for instance, by shaped or stochastic classical pulses. 

To unveil the true role of entanglement in virtual-state spectroscopy, we make use 
of two ingredients. First, we apply a full quantum formalism to the two-photon state, 
so we can identify clearly the amount of entanglement existing between the photons. 
Secondly, we consider a general form of the two-photon state, which allows us to consider 
different types of correlations and spectral shapes of the photons. 
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We will show that information about the energy levels of a medium by means of 
virtual-state spectroscopy can be obtained only when frequency correlations between 
the photons are present. Moreover, we will show that the correlations between them 
are not the only requisite for retrieving spectroscopic information of the medium. It 
is the combination of both, correlations and a specific spectral shape of the photons, 
which makes the realization of virtual-state spectroscopy possible. This result is of great 
interest because it specifies the type of two-photon source that needs to be used in order 
to experimentally perform the two-photon virtual-state spectroscopy technique. 



2. Light-matter interaction 

Let us consider the interaction of a medium with a two-photon optical field described 
by the interaction Hamiltonian Hi (t) = d(t) E^ (t), where d(t) is the dipole- moment 
operator and (t) is the positive-frequency part of the electric-field operator, which 
reads as E^ (t) = e[ +) (t) + E ( 2 +) (t). The electric field operators e{ +) (t) and E ( 2 +) (t) 
can be written as 



(t) = j 4^4« M ex P (-««>;*) - (!) 

where c is the speed of light, eo is the vacuum permittivity, A is the effective area of the 
field, and a (uj) is the annihilation operator of a photonic mode with frequency Uj. 

The medium is initially in its ground state \g) (with energy e g ). The probability 
that the medium is excited to the final state |/) (with energy Ef), through a two-photon 
absorption process, is given by second-order time-dependent perturbation theory as [IS] 

roo rti 



n 2 



/oo rt2 " 

dt 2 / dtMditi^Mzituh) , (2) 
-oo J — oo 



with 



M d (ti,t 2 ) = (fld^d^lg), (3) 

M^(t 1; t 2 ) = (%|£«(t 2 ) I*), (4) 

where |^//) denotes the final state of the optical field, which is assumed to be the vacuum 
state. 

Equation (|S]) can be expanded in terms of virtual-state transitions, to obtain 
Mi (h,h) = X; D (j) exp \-i (ej - i Kj /2 - E f ) t 2 ] 

x exp[—i(£g — £j + iKj/2)ti], (5) 

where = (f\ d\j) (j \ d\g) are the transition matrix elements of the dipole-moment 
operator. Equation (jSJ) shows that the excitation of the medium occurs through 
intermediate states \j), with complex energy eigenvalues £j — iKj/2, where Kj takes 
into account the natural linewidth of the intermediate states [19]. Also, we can write 
Eq. (0$ as 

(h, t 2 ) = E { 2 +) (t 2 ) E{ +) (*0 |^) + E[ +) (t 2 ) E { 2 +) (t x ) |^> , (6) 
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where we have kept the terms in which only one photon from each field contributes to 
the overall two-photon excitation. The first term of Eq. ([6]) corresponds to the case in 
which the photon field e[ + ^ (t) interacts first, and (t) interacts later. The remaining 
term describes the complementary case. 

Since we are interested in a two-photon absorption process, we consider the initial 
state of the optical field as an arbitrary two-photon state, which can be written as [20J 

|tf ) = J dv s dv& (u s , Ui) a\ (v s + ufj a\ (i/i + uf) |0) , (7) 

where s and i stand for signal and idler photonic modes, Uj = Uj — (j = s, i) are the 
frequency deviations from the central frequencies u®, and $ (u 8 , Vj) is the joint spectral 
amplitude, or mode function, which fully describes the correlations and bandwidth of 
the two-photon state. 

Finally, to quantify the degree of entanglement between the absorbed photons, we 
make use of the entropy of entanglement, defined as [21J 

E = -^A l log 2 A i , (8) 

% 

where \ are the eigenvalues of the Schmidt decomposition of the joint spectral 
amplitude, i.e., <&(v s ,Vi) = J2iV~Kfi(v s )g(i>i), with / and g corresponding to the 
Schmidt modes. It is worth remarking that the lack of entanglement between the pair 
of photons is characterized by a value of the entropy equal to zero. 

3. Two-photon absorption transition probability 

With the aim of recognizing in which situations virtual-state spectroscopy can be 
performed, we will compute the TPA transition probability of atomic hydrogen using 
different types of initial two-photon states. We have selected atomic hydrogen as model 
system, because it has been used in previous studies of virtual-state spectroscopy [18j [7] 
and it has been the subject of several one- and two-photon absorption experiments 
(221 |23l EH EE] . In our calculations, we will focus on the Is — > 2s two-photon transition. 
Due to quantum number selection rules [21] , this transition takes place via intermediate 
p states: Is —> {2p,3p, ...,np} — > 2s, which are coupled to the s states by real-valued 
transition matrix elements. The hydrogen atom energy levels are e n = — 13.6/n 2 eV 
(n = 1,2,3,...) and the natural linewidths of intermediate states Kj are taken from 
Refs. [221 El]- We assume the condition ej — e g = cu® + tuf, and that the final state 2s is 
Lorentzian broadened with a radiative lifetime = 122 ms [23], which is introduced 
in the model by averaging the TPA transition probability over a Lorentzian function of 
width K f [19]. 

3.1. TPA transition probability with uncorrelated classical pulses 

Let us consider first the case that has been studied in Ref. [T7]. It corresponds to 
the situation in which the two absorbed photons are embedded into rectangular-shaped 
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pulses of the same duration T p , with a tunable time delay r between them. This initial 
optical field can be represented by an uncorrelated two-photon state described by the 
normalized mode function 



T 

-^-sinc (T p u s /2) sine (T p Ui/2) exp [i (u a - v { ) r/2] . 

Z7l 



(9) 



With the state given by Eq. (j9]), and making use of Eqs. (j2J), © and dSJ), we can write 
the TPA transition probability as 

2 



<?->/ 



{Tp, r) 



where 



h 



h 2 e 2 c 2 A 2 T 2 

sin [Aw (T p - r) /2] sin [A f (T p - r) /2] exp [i A g (T p + r) /2] 



(10) 



AgAW 



A„A 



9^/ 



2i sin (A g T p /2) sin (A/t/2) exp [-i (A/T p - A ff r) /2] 



A g A / 



sin [Aw (T p - r) /2] sin [A f (T p - r) j2\ exp [zA g (T p - r) /2] 



A g Aw 



A„A, 



, (12) 



with Af = Ej — iKj/2 — Sf + u>q, A g = e g — Sj + iKj/2 + uo and Auj = e g — Sf + 2uo. For 
the sake of simplicity, we have assumed the condition w° = w° = co . 

We have computed the TPA transition probability for different values of T p and r. 
In all cases, it turns out be constant as a function of the delay (r) between the pulses, 
when t < T p , which implies that a Fourier analysis with respect to r would result in 
only one peak centered at zero-frequency, meaning that spectroscopic information about 
intermediate levels of the medium is not present in the TPA signal. 

From these results one can infer that when frequency correlations between the 
photons are not present, spectroscopic information about energy levels is not available, 
which implies that virtual-state spectroscopy cannot be performed by means of two 
delayed rectangular-shaped classical pulses, which is in contradiction with the results 
presented in [17] . 



3.2. TPA transition probability of a frequency correlated two-photon state with 
Gaussian spectral shape 

In view of the previous results, and the ideas and calculations presented originally in [7J, 
it naturally arises the question if the presence of frequency entanglement between the 
photons is the key ingredient that allows to access information about the energy level 
structure of the medium by means of two-photon virtual-state spectroscopy. In what 
follows, we will show that this is indeed not the case. 

In general, a two-photon state with tunable frequency correlations, and 
consequently tunable degree of entanglement, can be generated by means of type-II 
Spontaneous Parametric Down Conversion (SPDC), where two photons with orthogonal 
polarizations are generated in a second-order nonlinear crystal of length L, when pumped 






Figure 1. Joint spectrum of the two-photon state for different values of T_: (a) 
T_ = 2 ps, (b) T_ = 20 ps, and (c) T_ = 200 ps. (d) Entropy of entanglement as a 
function of T_ for Gaussian (red dashed line) and sine cardinal (blue solid line) shapes 
of the mode function. In all cases, the pump pulse duration is: T + — 10 ps. 



by a Gaussian pulse with temporal duration T+. After the crystal, signal and idler 
photons interchange their polarization and traverse a similar crystal of length L/2. After 
the addition of a tunable delay r between the photons, and restricting their spectrum 
using a Gaussian filter, the normalized mode function reads as 

/ rp rp \ 1/2 

-T 2 _(u s ~u t ) 2 /A 

X exp [iLNp (v s + Vi) /2 + zV;t] , (13) 

where TL = (N s — iVj)L/2, Nj (j = i,s,p) are the inverse group velocities. We have 
made use of the group velocity matching condition N p = (JVi + iV s ) /2 |26j, which ease 
the tuning of the frequency correlations, and the degree of entanglement, between the 
photons [27J. 

The frequency correlations of the down-converted photons can be tuned by carefully 
selecting the values of T + and T_. Figures 1(a) to 1(c) show the joint probability 
distribution of the two-photon state S (v s ,i>i) = |$ (z/ s , z/j)| 2 , which measures the 
probability of detecting a signal photon of frequency + v s in coincidence with an idler 
photon of frequency + U{. Frequency anti-correlated photons [Fig. 1(a)] [y s ~ —v^ 
are obtained when T + ^> T_, whereas for T + T_, we obtain frequency correlated 



± L _L 



exp 



-T, {i>i + v s 



exp 
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photons [Fig. 1(c)]. In the particular case when T_ = 2T + , frequency uncorrelated 
pairs of photons [Fig. 1(b)] are generated. Figure 1(d) (red dashed line) shows the 
dependence of the entropy of entanglement with XI for a fixed value of T+, for a mode 
function of the form given by Eq. ( TT3|) . 

Using the initial two-photon state described by Eq. (I13jl . we find that the TPA 
transition probability is given by 



P g ^ f (T_,T + ;T) 



32tTUJq r 2 . . 2 

X + X_exp -2T + {e g - E f + uj p ) 



h 2 e 2 c 2 A 2 



exp 



+F_ 



exp 



where r/^ = A^^ — iKj/2, with the energy mismatch given by 
the function F defined as 

2% K±- 



1 



7T JO 



exp (V) dy 



(14) 



e» — £„ — w , and 



(15) 



We have computed the TPA transition probability as a function of the delay 
between photons considering states bearing different types of correlations, particularly 
for uncorrelated and anti-correlated pairs of photons. As previously obtained, in the case 
of uncorrelated photons, the TPA signal is constant with the delay r, so no spectroscopic 
information about intermediate levels is available. 

Surprisingly, in the case of anti-correlated photons, the TPA transition probability 
is also constant with the delay r, which means that information about the energy level 
structure of the medium cannot be retrieved from the TPA signal either. This result is of 
great interest since it tells us that the use of a source of paired photons with entanglement 
alone is not enough to guarantee the successful retrieval of such information. Even 
though there exists entanglement between the photons that interact with the medium, 
there is another property of the two-photon state that is needed in order to perform 
virtual-state spectroscopy, namely a specific spectral shape of the frequency correlations. 



3.3. TPA transition probability of a frequency correlated two-photon state with a Sine 
cardinal spectral shape 

Fortunately, two-photon states with a Gaussian shape, which require a strong filtering 
of the pair of photons [28], are not naturally harvested in SPDC. By considering a 
more realistic shape of the mode function, we will show that two-photon virtual-state 
spectroscopy can retrieve the sought-after information about the energy level structure 
under a great variety of circumstances. 

As in the previous subsection, we consider a type-II SPDC process where 
an additional nonlinear crystal of length L/2 is used to achieve group velocity 
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Figure 2. Transition probability as a function of the delay r for anti-correlated 
photons (T_ = 2 ps). Pump pulse duration: T + = 10 ps. Y-axis in logarithmic scale. 



compensation. By introducing a tunable delay r between the photons, without 
restricting their spectrum, the normalized mode function is written as 

$(v a ,Ui)= - — 7= I exp -Tj (z/j + z/ s ) 2 sine [T_ - u t ) /2] 

x exp [zLA^p (z/ s + z/j) /2 + zz^t] . (16) 

The entropy of entanglement of the two-photon state described by Eq. (|T6l) is shown 
in Fig. 1(d) (blue solid line). Notice that in this case, due to the presence of the sine 
cardinal function, only quasi-uncorrelated photons can be generated. 

We now make use of the initial two-photon state described by the mode function 



given in Eq. (1161) to write the TPA transition probability as 

'V2T. 



P g ^ f (T_,T + ;r) 



647ra;g 



h 2 e 2 c 2 A 2 T„ 



71 



exp 



2T^(e g -e f + u Jp y 



exp 



-iri® (T. 



exp 



-i7] ij) (T_ + r)] } 



(17) 



where = D^/rj^. 

Figure 2 shows the TPA transition probability as a function of the delay between 
the pulses. Notice the nonmonotonic behavior of the TPA transition probability for 
anti-correlated photons. This means that spectroscopic information is contained within 
the TPA signal, which might be related to the energy level structure of the medium. In 
order to retrieve this information, we follow [7] and perform an average of Eq. (II 7p over 
a range of values of T_ to obtain the weighted-and-averaged TPA transition probability 

P(r) = - j Tm . m P g ^ f (T_, T + ; r) T_dT_, (18) 

where T = T max - T min . 



4 5 6 7 8 9 

A< 3) (eV) 



Figure 3. Fourier transform of the normalized weighted- and- averaged TPA transition 
probability as a function of the energy mismatch A"'. The delay range considered is 
< t < 2 ps, with an integration time of 2 < T_ < 10 ps. Y-axis in logarithmic scale. 



Provided that T 1/ A^') — , to eliminate unwanted terms at intermediate 
frequencies, a straightforward Fourier analysis of Eq. ffl8|) reveals the curve shown in 
Fig. 3. We see that peaks emerge from the Fourier transform of the weighted-and- 
averaged TPA transition probability, whose locations determine the energy mismatch 
of the intermediate states: 5.1, 6.98, 7.65, 7.95, and 8.12 eV. With these values and 
the definition of the energy mismatch, we obtain the virtual-state energy values: —3.40, 
— 1.51, —0.85, —0.54, and —0.37 eV. These energy values can be readily identified with 
n = 2, 3, 4, 5, 6 corresponding to the 2p, 3p, 4p, 5p, and 6p states, respectively. 

To obtain a clearer image that two-photon virtual-state spectroscopy depends on 
the quantum interference from different contributions of intermediate-state transitions 
with a specific spectral shape, let us consider a simpler, even though ideal, case where 
a single intermediate quantum state (3p, Ap or bp) is present [28]. Figure 4 shows the 
two-photon transition probability as a function of T_ for a fixed value of T + and delay 
t = 0, considering three different intermediate states. Changing T_ modifies the degree 
of entanglement between the photons (as shown in Fig. 1). Figure 4(a) corresponds 
to the case of frequency entangled photons with a gaussian spectral shape, while Fig. 
4(b) depicts the case of a sine cardinal spectral shape. From these results we find that, 
even though in both cases the degree of entanglement between photons is the same, 
when considering frequency entangled photons bearing a sine cardinal spectral shape, a 
nonmonotonic behavior of the TPA transition probability appears, whereas in the case 
of a Gaussian spectral shape this is no longer true. 

By making use of entangled photons bearing a sine cardinal spectral shape we can 
then retrieve information about the atomic transitions by means of a Fourier analysis of 
the weighted-and-averaged TPA signal [Eq. ([TBI ]. It is worth mentioning that the same 
information as the one depicted in Fig. 3 can also be obtained when quasi-uncorrelated 
photons (T_ = 40 ps) are used, meaning that virtual-state spectroscopy can also be 
performed even with a low degree of entanglement between the photons. 
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Figure 4. Single-intermediate level TPA transition probability as a function of T_ for 
(a) Gaussian mode function and (b) sine cardinal mode function. Intermediate levels 
correspond to: 3p (solid line), 4p (dashed line) and 5p (dotted line). Delay is set to 
t = ps. Pump pulse duration: T + = 10 ps. Y-axis in logarithmic scale. 



4. Conclusions 

We have shown that two-photon virtual-state spectroscopy cannot be performed by 
means of two uncorrelated rectangular-shaped classical pulses. This implies that, in 
order to extract information about the energy levels of a medium, one has to make 
use of light bearing frequency correlations. Notwithstanding, we have found that these 
are not the only requisite for accessing such information. In fact, it is a combination 
of two-photon frequency correlations and a specific spectral shape of the correlations 
which allows one to perform two-photon virtual-state spectroscopy. 

The results presented here are of great interest because they help to end the 
discussion whether entanglement is needed or not, and also, to identify clearly which 
types of two-photon sources can be used to experimentally implement virtual-state 
spectroscopy. By clarifying the role of entanglement, we have found that even though 
it can be required, not any source of entangled paired photons can be used to 
experimentally perform such technique. 
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